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ABSTRACT 


Natural convection liquid cooling of simulated electronic compo- 
nents in a vertical channel was investigated. The test surface 
contained a single column of eight rectangular, protruding heated ele- 
ments, each simulating a 20-pin dual-in-line package. Temperature 
measurements and flow visualization were performed for a number of 
power dissipation levels and channel widths. Collectively, this infor- 
mation was used in interpreting the flow and transport characteristics. 
A correlation to predict the heat transfer rates was developed based 
on the component surface temperatures. Optimum channel widths 
were determined from these surface temperature measurements for 
the range of power levels investigated. Temperature distributions in 


the fluid were measured using a traversing thermocouple probe. 
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I. INTRODUCTION 


A. STATEMENT OF PROBLEM 

The continuing trend toward microminiaturization of electronic 
components has made thermal design critical to achieving higher 
packaging densities. The ever-increasing heat fluxes must be removed 
while keeping the device junction temperatures typically below 100°C 
[Ref. 1]. The failure of microelectronic chips increases exponentially 
as the junction temperatures increase [Ref. 2]. 

Microchips act as discrete sources of thermal energy. Effective 
removal of the generated energy is of major import in determining 
how they might be used without overheating the elements themselves 
or the surrounding package. Thermal control is one of the critical 
considerations in the design of any electronic cooling system. The 
thermal designer’s task is further complicated by the presence of 
arrays of electronic components placed in proximity to each other, 
requiring detailed considerations of transport interactions between 
neighboring elements. 

The design of packaged electronic components must consider all 
anticipated modes of heat transfer. A complete thermal analysis will 
account for all three mechanisms of thermal energy transport: 

¢ conduction 
¢ convection 


e radiation 


The areas of concern include transferring heat, via conduction, from 
internal heat sources to the external surface of the package and the 
subsequent removal of heat from the external surface by a cooling 
medium. Also, the coolant temperature must be maintained at pre- 
scribed levels for a given environment. 

The present work focuses on the heat removal from electronic 
packages through convective cooling. This convection may, in general, 
be forced or natural; the medium liquid or gas. Because of the higher 
density and thermal conductivity of liquids compared to air, heat 
transfer coefficients an order of magnitude or more larger than for air 
are attainable. Due to the significantly high heat transfer rates achiev- 
able, cooling of electronic components by direct liquid immersion has 
been identified as an extremely effective thermal control technique 
[Ref. 2]. 

Convection heat transfer rate from a surface at temperature T; of 
area A to surrounding fluid environment Tyyf is given by Newton's law 


of cooling: 
q = AN(T, - Tee) 


This is a definition of h, the convective heat transfer coefficient, more 
than a phenomenological law. Evaluation of h depends upon heat 
source geometry, orientation, fluid properties, and fluid velocity. 
Natural convection results when fluid density gradients exist in 
the presence of a body force field such as gravity. Density gradients in 


the fluid commonly occur due to temperature gradients. Advantages of 


employing natural convection in electronic cooling include high relia- 
bility, low cost, reduced noise, and a simpler design. Despite these 
potential benefits, only a limited number of studies are available on 
liquid immersion natural convection cooling of electronic components. 

Most existing natural convection studies, both experimental and 
analytical, have considered heating surfaces that were semi-infinite. 
Electronic component cooling configurations typically employ arrays of 
small, discrete, flush, or protruding heat sources. The available heat 
transfer correlations are inadequate in predicting transport in these 
geometries. Convective heat transfer from an array of microchips is 
strongly dependent on package shape and array density. Indeed, in 
real systems, individual elements or modules often vary significantly in 
size, shape, and in the manner in which they are arranged on the 
substrate. 

The geometric complexity is further compounded by the nonuni- 
form nature of heat dissipation from the array. Individual component 
dissipation is a function of its utilization. Variation in power dissipa- 
tion from component to component on the same board and from board 
to board within the system usually exists. The orientation of the 
boards on which the elements are mounted is yet another variable. 

Often, these semi-regular arrays face the back of an adjacent 
board. The resulting vertical channel has one surface which is rela- 
tively smooth and the other one covered by an arrangement of large, 
heat-dissipating elements. Within the channel, local convective heat 


transfer is driven by two phenomena, a local buoyancy force due to the 


heated components and a forced convection effect arising from the 
flow history, over the entire channel length, including both heated and 
unheated regions. A simultaneous consideration of these mechanisms 
is essential in obtaining predictive capability for transport rates in 
these systems. 

The complexity of the problem described demands an approach 
which builds on the fundamental physical processes. Efficient cooling 
cannot be attained without understanding the heat transfer from each 
package by a determination of its flow and thermal fields. The overall 
goal is to develop a methodology based on sound physical principles 
which will allow the prediction of individual component temperatures, 


by superposition, in a nonuniformly heated array of components. 


B. IMMERSION COOLING—ANALYTICAL AND EXPERIMENTAL 
STUDIES 


Reference 3 provides an excellent literature review of previous 
work, both analytical and experimental, in the field of liquid immer- 
sion cooling of electronic equipment. Additionally, very relevant work 
has been performed by Ortega and Moffat [Refs. 4-6]. They have 
experimentally examined natural convection air cooling of simulated 
electronic components protruding from a vertical surface, with and 
without a shrouding wall. Heat transfer rates significantly higher than 
for a smooth vertical plate were found. They interpreted this as an 
indication of turbulent free convection. When a shroud was in place, 
they found plate-averaged heat transfer coefficients 40 to 50 percent 


higher than for a smooth, parallel plane channel. 


Jaluria [Ref. 7] performed two-dimensional computations for 
natural convection air flow over a vertical surface with multiple wide, 
flush-mounted heaters. The study treated the flow as a boundary layer 
problem. Rajakumar and Johnson [Ref. 8] reported computations of 
free convection heat transfer using rectangular strip heating surfaces. 
Their finite element approach showed significant stagnant areas 


between the blocks. 


C. OBJECTIVE 

An experimental study was undertaken to to determine the 
hydrodynamic and thermal characteristics of buoyancy-induced flow 
over several heated protruding components mounted on a vertical 
surface. 

A single column array was chosen because it provided the sim- 
plest geometry featuring discretized heat dissipation. The effects of a 
neighboring card were simulated by a movable vertical shroud plate. 
Specific objectives of this study were: 

e to visualize the transient and steady-state, natural convection flow 
within an interrupted channel, for a range of component power 
dissipation levels and channel widths; 

e to measure component temperatures for various power inputs and 
shroud spacings to develop an appropriate nondimensional cor- 
relation to predict the heat transfer rates for this geometry; and 

e to measure temperature distributions in the adjacent fluid, using a 


traversing thermocouple probe to assist in evaluating the hydro- 
dynamic and thermal nature of the flow. 


This study was a follow-on to work described in Reference 3 and 
is part of a continuing effort to study natural convection liquid immer- 


sion cooling of electronic equipment. 


Il. EXPERIMENT 


A. GENERAL DESIGN CONSIDERATIONS 

Detailed descriptions of equipment design and construction are 
available in Reference 3. Only a brief overview is presented here. 

As seen in Figure 1, the assembly consisted of a vertical test sur- 
face which contained a single column of eight rectangular protruding 
heated elements. Each protrusion simulated a 20-pin dual-in-line 
package. The elements were stainless steel blocks spaced one inch 
apart, center to center, on a plexiglass substrate. A 12.70 mm thick 
plexiglass shroud, when used, simulated the reverse side of a second 
printed circuit board. 

The blocks (Figures 2 and 3) had imbedded thermocouples on 
each of the five surfaces in contact with the fluid to provide surface 
temperature measurements. A sixth thermocouple was placed in the 
center of each element mounting slot for measurements of heater 
temperature. A miniature thermocouple probe mounted on a three- 
dimensional traverse enabled measurements of temperature distribu- 
tions within the fluid. 

A foil heater mounted on the component side facing the substrate 
was powered using a regulated D.C. power supply. Each heater was 
run in series with a precision resistor and all eight heaters were in 


parallel with the power supply. The current to each heater was 
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calculated by subtracting the heater voltage from the source voltage, 
then dividing by the resistance of the precision resistor. Power dis- 
sipation through each element was determined by calculating the 
product of heater voltage and heater current. 

Figure 2 shows the orientation of the coordinate axes used in the 
discussion to follow. The front face of the component is that surface 
in the x-z plane in contact with the fluid. The side faces are the two 
x-y surfaces, right and left as determined by viewing the component 
from the front. The top and bottom faces are those in the y-z plane, 
downstream and upstream of buoyant fluid flow, respectively. 
Component height is in the x-direction, depth in y, and breadth in the 
z-direction. 

The immersion fluid was purified water contained in a one cubic 
meter, plate glass tank. Suspended in the water were particles of 
Pliolite, an inert pigment with a specific gravity of 0.93. The particles 
were illuminated with an eight milliwatt helium neon laser split into a 
plane by use of a cylindrical lens. The plane of visualization was varied 


by realignment of the laser lens assembly and the camera. 


B. EXPERIMENTAL PROCEDURE 
1. Established Method 
Details are again provided in Reference 3, hereafter referred 
to as the previous study. Noteworthy differences are noted here. 
2. Instrument Settings 
Power settings of 0.2, 0.5, 1.0, 1.2, and 1.5 watts were used 


for each channel width. Because of the inability to de-aerate the 
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immersion fluid, 1.5 watts was the upper limit of power available. 
Settings above this generated bubbles on the block surfaces. 
3. Instrument Readings 

When measuring steady-state surface temperatures, succes- 
sive temperature measurements at 10-minute intervals were com- 
pared. Only when all thermocouples varied less than 0.10° C were the 
final measurements taken. The temperature acquisition program of 
the previous study, modified to reflect the numbering of elements 
from bottom to top (flow direction) and to allow data storage 
(Appendix A), was then executed. 

When using the traversing thermocouple probe, the data 
acquisition system was programmed to continuously monitor the volt- 
age registered by a single channel. Observation of the on-screen print- 
out allowed determination of temperature variation from point to point 
within the flow field, as well as at any point over time. 

A program (Appendix A) which monitored from one to eight 
thermocouples over time was used when making a record of the tran- 
sient temperature response of the surface thermocouples to power 
step up or step down. Typical data presented were for the front face 
temperatures of the various components. 

4. Photographic Technique 

Best results were obtained when using an f stop of 2.8 and 
exposure time of 20 to 40 seconds. Camera placement was approxi- 
mately six inches from the tank wall at a height which allowed all 


eight blocks in the field of view. The camera was inclined at an angle 


eZ 


of about 5 degrees from the test surface when a shroud was present. 
This was helpful in minimizing reflected light. A shroud spacing of 
two times component height from the plexiglass substrate was deter- 


mined to be the minimum at which useful pictures could be obtained. 


C. DATA ANALYSIS 

While Qconp and Qconv were calculated exactly as recommended 
by the previous study, further data reduction differed significantly. 

For each block, a surface-averaged temperature was generated by 
multiplying the temperature and area of each fluid-exposed face and 
dividing the summation over the block by the total exposed face area 
of the block. The average over the eight blocks of the surface-averaged 
temperatures was then averaged with the immersion fluid tempera- 
ture to produce a film temperature. Fluid properties based on curve 
fits (Appendix B) of data from Incropera and Dewitt [Ref. 9] were then 
calculated for the film temperature. 

A new characteristic length, L, for determining a nondimen- 


sional temperature excess, a a modified Grashof number, Gr*, and 


0 
Nusselt number, Nu, was chosen. The progression of its use for each 


block is as follows: 
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The various symbols used are defined in the List of Symbols section. 
Other characteristic lengths, L to determine block spacing along 
the plexiglass substrate and D to represent the horizontal breadth of 


the element, were also chosen. Their use is as follows: 
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a = nondimensional location of component centers along the 
test surface 
a = horizontal position along the element's front face 


Quantitative development was best expressed as plots of: 
e Nu vs. Gr* relationships 


e nondimensional block surface temperature excess vs. element 
position 


e temperature distributions in the fluid adjacent to the test surface 


e block surface temperatures vs. time 


Pertinent trends in these variations are diScussed next. 
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I. NO SHROUD—INFINITE CHANNEL WIDTH 


A. ALL ELEMENTS POWERED, STEADY STATE 
1. Flow Visualization 

Photographs of the natural convection flow in a plane (x-y) 
geometrically centered through each component are shown in Fig- 
ure 4. The visualizations represent all elements heated to power 
settings of 0.2, 0.5, and 1.0 watts, left to right, respectively. The 
exposure times were 40, 20, and 20 seconds, left to right. The pho- 
tograph toning techniques required for mass printing of this work do 
not allow sufficient resolution for detailed study. Should the reader 
desire to see prints of much greater clarity, contact Y. Joshi at the 
address given in the distribution list of this study. 

The flow appears to consist, in each case, of strong upward 
motion driven by the buoyant force near the test surface. The fluid 
tends to follow the contours of the protrusions. The particle traces 
are closer passing the element front face as compared to the inter- 
block space. Dead zones are not obvious except perhaps for very small 
areas where the blocks’ top faces intersects the substrate and for a 
slightly larger area above the top face of the last block in the array. 
There is no evidence of recirculation or formation of dead zones near 


the bottom face of any block. 
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Figure 4. Steady Flow in X-Y Plane at Power Levels of 0.2 


and 1.0 Watts. Powers Increase From Left to Right 
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In all three photographs, a developing boundary layer characteristic is 
evident farther out from the protrusions. It attains a thickness of 
about 3.5 component depths for 0.2 watts and 2 for 1.0 watt. Its 
origin moves upstream with increased power, from about 1 com- 
ponent height at 0.2 watts to 3 component heights at 1.0 watt. The 
velocities appear greater and the boundary layer becomes thinner as 
power is increased. An. ent fluid is entrained over the entire length 
of the layer. 

Figure 5 presents flow visualization in the plane (x-z) passing 
through the front faces of the elements. All elements were powered 
to 1.2 watts and exposure time was 40 seconds. The photograph 
shows an upward-moving buoyant layer near the components. Ambient 
fluid is entrained from each side into the flow, resulting in an increase 
in the horizontal extent of the upflow region downstream. The flow in 
the vicinity of the elements appears three-dimensional. 

As with the side view, the flow appears to follow the contours 
of the projections to about the same degree (dipping in about 3-5 mm 
from the side faces). the flow also appears to wrap around the cor- 
ners. The particle traces about 20 percent along the block (z direc- 
tion) from the side edge appear to be largely in the x-y plane and 
those at the side face are in the x-z plane. The particle traces in 
between are 3-D, each progressively moving along the 90° arc. 

Overall, the flow can best be described as following the con- 


tour of an hourglass split in the x-z plane. The elements are the upper 
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gure 5. Flow in X-Z Plane, 1.2 Watts 
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and lower vessels. This configuration is repeated over the length of 
the array. It becomes more pronounced as the flow moves up the 
array and is eventually shed as a plume above the last element. 

2. Quantitative 

Figures 6 and 7 are nondimensional representations of tem- 
perature excess for 0.2 and 1.0 watts, measured on various component 
surfaces vs. block position. Also shown are the data for heater 
temperature measurements. The two plots follow the same patterns. 
The numbers for the 0.2 watts are, as expected, larger due to the 
inverse relationship with T*, which is h dependent. Also as antici- 
pated, the heater excesses are much greater. The drop in the heater 
temperature between blocks 6 and 7 could perhaps result from a 
slightly smaller (~4% less) power achieved by block 7. Unfortunately, 
the heater thermocouple for block 8 was inoperative. For both power 
settings, the imbedded surface thermocouples showed increasingly 
greater excesses along the array. 

Of primary importance in thermal design is the availability of 
the appropriate heat transfer relationships. As discussed earlier, such 
correlations for discrete protruding heat sources in liquids have not 
been available in literature up to now. To this end, a logarithmic plot 
of Nu vs. Gr* (Figure 8) was generated with eight curves representing 
element position. Each curve had five points for the power settings 
0.2, 0.5, 1.0, 1.2, and 1.5 watts resulting in a range of Gr*. Although a 


linear pattern with Gr* was evident, the spread in data indicated an 
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Figure 6. Nondimensional Temperature Excess vs. 
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Nondimensional Temperature Excess vs. 


Figure 7. 
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Figure 8. Data for Nusselt Number vs 


additional dependence on element location. A review of the data sug- 
gested a method by which it could be collapsed. A modified Nusselt 
number, Nu*, was calculated by use of a scaling factor based on the net 
convected energy of the flow at a downstream location. Since each 
block dissipates approximately an equal amount of power, the net 
convected energy is proportional to the number of upstream compo- 
nents. The scaling factor was of the form Bn@ where Bn is the block 
number and a>O. A value of a= 1/6 was found to provide the best 
correlation for the data. As expected, the strong three-dimensional 
effect on transport results in a weaker downstream dependence than 
the 1/5 exponent in the Fujii and Fujii correlation for a uniform flux, 


semi-infinite vertical surface. Therefore, 


Nu*= Nu- Bn‘ 


for each element at each power setting. The modified data are plotted 
in Figure 9 along with the linear least squares fit. Based on this, the 


following correlating equation for the heat transfer rates is obtained: 


= be 
Nu*= 188(Gr *) 


for 50 < Gr, < 1000 and Pr ~ 7. 
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B. ALL ELEMENTS POWERED, TRANSIENT 
1. Flow Visualization 

Figure 10 is a series of photographs taken during the 0.2 watt 
step-up transient. For each picture, the exposure time was 20 sec- 
onds. Exposure of the first picture began with power application. The 
others were initiated at 20, 80, 120, and 200 seconds, respectively, 
following the start of heating. Clearly, the first two show predom- 
inantly parallel flow. The next pictures represent increasing devel- 
opment of the contour flow and boundary layer characteristic 
development. By 200 seconds, almost fully developed flow is 
recognized. 

Figure 11 is a similar series, again of 20-second exposures, 
taken during the step-down transient from 1.0 watt. The first visual- 
ization at far left began at power down. The remaining photographs 
were initiated 80, 160, 240, and 440 seconds, respectively, following 
shutdown. The steady-state flow characteristics persist to an 
increasingly lesser degree for the entire period observed. The 
entrainment of ambient fluid lessens at large times. Additionally, flow 
origin moves slightly downstream over time. 

2. Quantitative 

The 1.0-watt step-up, step-down transient temperature 
responses for the component front surfaces are presented on the 
same graph (Figure 12). Steady state is achieved about 25 percent 
faster for the power up. The curves representing blocks 1 through 8 
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Figure 10. 0.2 Watt Step-Up Transient During 
Various Time Intervals 
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1.0 Watt Step-Down Transient Durin 
Various Time Intervals 


Figure 11. 
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Figure 12. Transient Temperature Response Following 
Initiation of Heating, 1.0 Watt 


JES) 


excluding 3 travel closely together, especially in the first quarter of 
the response period for the step. The almost single line corresponds 
to the parallel flow pattern seen in the photographs. This indicates a 
purely diffusive transport response for a large part of the transient. 
Later, convective effects become appreciable as entrainment develops 
downstream. For the step down, steady state temperature is reached 
much more quickly than breakdown of the established flow. The per- 
sistence of flow is conducive to the rapid cooling of the blocks. 

The plot for 0.2 watts (Figure 13) reaches steady state at 
about the same time for both step up and step down. This period is 
longer than for either 1.0-watt case and indicates that the stronger 


flow more readily establishes the final thermal pattern. 


C. TEMPERATURE MEASUREMENTS IN THE FLUID 
Characteristics of the thermal transport in the fluid were mea- 
sured using a traversing thermocouple probe adjacent to all protruding 
elements and for points midway between the components (Figures 14 
and 15). The fluid temperature responses near the blocks and inter- 
block spaces follow the same pattern, with temperature levels in the 
vicinity of blocks 4 and 6 surprisingly greater than for the others. 
Indeed, both graphs can be seen to have four groupings. Block 1 by 
itself, the lowest curve, is followed by 2, 3, and 5; then 4, 7, and 8; 
and finally block 6 by itself. More detailed local measurements are 


needed to explain these variations. 
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Figure 14. Fluid Temperature Distribution in the Normal 
Direction, Away From Block Faces 
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The layer between the blocks is slightly thicker with respect to 
that at the block faces. This might be expected from the flow 
visualization (Figure 4) showing that the particle traces spread out 
between elements. In each case, the thermal layer is about one 
component height or less and is therefore deeply imbedded within 
the momentum exchange region. This is expected due to the 
relatively large Pr ~ 7 for water. 

Next, the traversing probe was placed in the thermal layer, 1 mm 
from the surface, and moved in the x-z plane across the test surface 
face (Figures 16 and 17). The temperature variations, with position, 
reinforce the photographs in suggesting strongly three-dimensional 
transport. Ambient temperature is the lowest point shown. The 
peaks at about 15 percent from each edge (Figure 16) correlate with 
the regions in the visualization (Figure 5), where the flow begins to 
shift from x-y dominant to x-z and may indicate a common region 
through which those particles are flowing. The result is more, 
warmer fluid convected over that region. Likewise, the peaks of Fig- 
ure 17 might represent a similar flow shift at locations between the 
components. 

Finally, the probe was placed above the center of the top face of 
block 8, 2mm out from the substrate to measure the downstream 
temperature decay (Figure 18). The extremely high temperature in 


the immediate vicinity of the top face suggests a dead zone, as did the 
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Figure 18. Fluid Temperature Decay Downstream 


of Top Component 
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visualizations, where a small pool of fluid is heated primarily by ther- 
mal conduction from the component upper surface. There was no 
sign of a vortex. 

All measurements with the probe were remarkably consistent 
over time (+1 pv) at any position. No evidence by any measurement 
indicated the presence of turbulent flow for the power levels 


investigated. 


D. INDIVIDUAL ELEMENTS POWERED 

Individual blocks were powered to 1.0 watt, one at a time. The 
resulting steady-state array temperatures were monitored and the flow 
visualized. These are discussed next. 

1. Flow Visualization 

A series of photographs show single-element heating as that 

element changes from one block to the next up the array. The expo- 
sure times for all pictures presented were 40 seconds. Figure 19 pre- 
sents left to right blocks 1, 3, and 4 heated and Figure 20 shows 
blocks 5, 6, and 7. The flow with only the lowest block heated resem- 
bles the all-components-heated case, indicating flow around block 
contours and a boundary layer characteristic over the entire length. As 
progressively higher blocks are heated, a quiescent area begins to 
develop adjacent to the lower part of the test surface. Typically, fluid 
entrainment results in vigorous flow around the first block upstream of 
the one heated. At two blocks upstream, a comparatively weaker flow 


is seen and at three the fluid is nearly quiescent. 
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2. Quantitative 


Nondimensionalized temperature excesses for the faces and 
heater of each element are displayed (Figures 21 through 25). The 
plots show very similar temperature excess patterns around the 
heated element as the selected component progresses up the array. 
The protruding elements have minimal effect when they are part of 


the unheated starting length. 


E. COMBINATIONS OF ELEMENTS HEATED 

As in Section D, progression was made up the array, except that 
all upstream blocks remained heated. Figures 26 through 30 indicate 
very little effect of powering of downstream components on the tem- 
perature responses of already-powered upstream components. Also 


evident is a very sharp decay in temperature in the unheated blocks. 
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Figure 25. Nondimensional Temperature Excess Levels at 
Heaters of Individually Heated Components 
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Figure 28. Nondimensional Temperature Excess Levels at 


Top Faces of Multiple Heated Components 
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Figure 29. Nondimensional Temperature Excess Levels at 
Bottom Faces of Multiple Heated Components 
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Figure 30. Nondimensional Temperature Excess Levels at 


Heaters of Multiple Heated Components 
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IV. SHROUD IN POSITION—FINITE CHANNEL WIDTH 


A plexiglass plate (shroud) was used (Figure 1) to simulate the 
channel formed by two adjacent parallel printed circuit boards. The 
test surface as before has protruding heated elements. The shroud 


simulates the reverse side of the second board in an array. 


A. ALL ELEMENTS POWERED, STEADY STATE 
1. Flow Visualization 

Photographs presenting natural convection flow in the x-y 
plane, geometrically centered through each component (Figure 31), 
were evaluated. The visualizations are for all elements heated to 0.5, 
1.0, 1.2, and 1.5 watts. The exposure times were 20 seconds. The 
shroud distance was 18 mm from the substrate or approximately three 
component depths. The first six blocks are clearly visible in each pic- 
ture. Figure 32 shows the top portion of the test surface for 0.5 watts. 
The other power settings do not differ significantly and are not shown. 

The flow near the elements appears similar to that observed 
in the absence of the shroud. However, the shrouding wall prevents 
ambient fluid from entrainment into the buoyant up-flow in this plane. 
As power is increased, the paths described by the particle traces are 
more closely spaced, indicative of increased velocity and increased 


fluid flow. 
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Figure 32. Steady Flow in X-Y Plane, 0.5 Watts, 18 mm 
Channel Spacing, Top Blocks 
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Figure 33 is the front view of flow in the x-z plane through 
the front surface of the elements. It is for 1.0 watt, 18 mm shroud 
spacing, and exposure time of 40 seconds. Unlike the no-shroud con- 
dition, flow is more restricted to the sides of the protruding blocks. 
The characteristic hourglass flow seen for infinite channel width is 
weak in the vicinity of the middle two blocks and absent elsewhere. 
The entrained fluid turns very sharply into the flow direction. Because 
the flow is restricted to the channel width and is predominant on the 
block sides, its velocity is much greater here. Consequently, there is 
more forceful entrainment of fluid in the x-z plane. 

2. Quantitative 

Due to the large amount of data collected, temperature 
excesses for power settings of 1.2 and 0.2 watts only are presented in 
graphical form. These settings are representative of the upper and 
lower ends of the power range available. Channel widths presented 
are 6.2, 8, 14, 18, and 23 mm measured from the substrate. 

The first width was just barely greater than component 
height. This spacing is representative of pure conduction dominated 
transport. The nondimensionalized temperature excess vs. element 
position graphs for 0.2 watts (Figures 34-38) and 1.2 watts (Figures 
39-43) show significant difference between the smallest spacing and 
all others. With that spacing excluded, Figures 44 through 53 present 
the relevant data in a manner better suited for analysis. The plots for 


0.2 watts show the next closest spacing, 8 mm from the substrate, to 
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Figure 33. Flow in X-Z Plane, 1.0 Watt, 18 mm Channel Spacing 


96 


t 
L 
U 
U 
I 
T 
I 
I 
4 
L 
’ 
rT 
I 
J 
$ 
$ 
U 
T 
t 
‘ 
U 
$ 
$ 
r 


GNouHs O 
WW te 
WW vl 
WW VI 
WIN 8 
WW cg 
(NADA T 


© 
~ 


d+ 


O 
O 





weweewewenee we eweweeweeweewe nw w= wee = = 





ee ee ee 


008 OS OOF ose coe ogz 
e-J7 «LALd 


pcs 


2 Watts 


Figure 34. Nondimensional Temperature Excess, O 


Various Spacings, Front Faces 
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Figure 35. 
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Figure 36. Nondimensional Temperature Excess, 0 
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Figure 37. Nondimensional Temperature Excess, 0.2 Watts 
Various Spacings, Bottom Faces 


60 


| 
1 
' 
4 
$ 
Tr 
J 
$ 
$ 
$ 
$ 
T 
$ 
U 
' 
$ 
$ 
T 
J 
4 
U 
J 
| 
Tr 
t 
| 
4 
J 
' 
T 
U 
$ 
U 
$ 
' 
T 
J 
J 
i] 
' 
' 


w-ererereeeerereeewereewe sew eweewe wewew 2 wee = 


= 





2 Watts 


Figure 38. Nondimensional Temperature Excess, 0 


Various Spacings, Heaters 
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Figure 39. Nondimensional Temperature Excess, 1.2 Watts 
Various Spacings, Front Faces 
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Figure 40. Nondimensional Temperature Excess, 1.2 Watts 
Various Spacings, Right Faces 
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Nondimensional Temperature Excess, 1.2 Watts 
Various Spacings, Top Faces 


Figure 41. 
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Nondimensional Temperature Excess, 1.2 Watts 
Various Spacings, Bottom Face 


Figure 42. 
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Figure 43. Nondimensional Temperature Excess, 1,2 Watts 
Various Spacings, Heaters 
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Figure 44. Nondimensional Temperature Excess, 0 


Excl. 6.2 mm, Front Faces 
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Figure 45. Nondimensional Temperature Excess, 0.2 Watts 
Excl. 6.2 mm, Right Faces 
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Figure 46. Nondimensional Temperature Excess, 0 
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Nondimensional Temperature Excess, 0.2 Watts 
Excl. 6.2 mm, Bottom Faces 


Figure 47. 
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Figure 48. Nondimensional Temperature Excess, 0 


Excl. 6.2 mm, Heaters 
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Figure 49. Nondimensional Temperature Excess, 1.2 Watts 
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Figure 50. Nondimensional Temperature Excess, 1.2 Watts 
Excl 
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Nondimensional Temperature Excess, 1.2 Watts, 


Figure 51. 


Excl. 6.2 mm, Top Faces 
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Nondimensional Temperature Excess, 1.2 Watts, 


Figure 52. 


Excl. 6.2 mm, Bottom Faces 
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Figure 53. Nondimensional Temperature Excess, 1.2 Watts 
Excl. 6.2 mm, Heaters 
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yield higher temperatures over the range of element positions. For 
higher spacings, there was no consistency in the relative position of 
the curves from one element position to the next. This indicates 
relative independence of transport on channel width for spacings 
greater than 8mm. However, data for 1.2 watts includes even the 
8 mm spacing within the mix of curves (its curve does not consistently 
yield higher temperatures). Therefore, for 1.2 watts the indepen- 
dence of transport extends below the 8 mm spacing. General patterns 
observed for both powers are: 


e a steady rise in front and side face temperatures with element 
position as expected due to exposure to increasing convected 
energy at higher positions. 

e a relative levelling of top and bottom temperatures at higher posi- 
tions. This is perhaps due to increasing fluid velocities along the 
channel. 


e the drop in heater temperatures at block 7, which has four per- 
cent less power on average. 


The differences between the two powers result from the 
stronger driving force at 1.2 watts. Increased fluid velocities result in 
improved heat transfer characteristics, making 8 mm an acceptable 
spacing. This point is emphasized by referring to the graphs of com- 
ponent face temperature excess vs. shroud spacing for 0.2 watts 
(Figures 54-58) and 1.2 watts (Figures 59-63). Significantly, for the 
1.2 watts there is no further temperature reduction for a spacing 
greater than 8mm. _ Indeed, there is a slight rise in temperature of 


the latter blocks (7 and 8) when the shroud is removed. These higher 
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Figure 54. Nondimensional Tem erature Excess vs. Shrou 


Position, 0.2 Watts, Front Faces 
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Figure 55. Nondimensional Temperature Excess vs. Shroud 


Position, 0.2 Watts. Right Faces 
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Figure 56. Nondimensional Temperature Excess vs. Shroud 
Position, 0.2 Watts, Top Faces 
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Nondimensional Temperature Excess vs. Shroud 


Figure 57, 
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Figure 58. Nondimensional Temperature Excess vs. Shrou 
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Figure 59. Nondimensional Temperature Excess vs. Shroud 


Position, 1.2 Watts, Front Faces 
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Figure 60. Nondimensional Temperature Excess vs. 


Position, 1.2 Watts, Right Faces 
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Figure 61. Nondimensional Temperature Excess vs. Shroud 
Position, 1.2 Watts, Top Faces 
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Nondimensional Temperature Excess vs. Shrou 
Position, 1.2 Watts, Bottom Faces 


Figure 62. 
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Figure 63. Nondimensional Temperature Excess vs. Shrou 
Position, 1.2 Watts, Heaters 


ow 


temperatures (at the last blocks) are not observed when a shroud is in 
position because the global forced convection within the channel is 
driving the fluid more strongly past the last block. As seen from Fig- 
ures 54 through 58, face temperatures at a power level of 0.2 watts for 
8 mm are somewhat higher than for greater spacings. By 14 mm, the 
difference with larger spacings begins to become indistinguishable. 
Clearly for this geometry and for the power settings examined, a 
channel width of 2.5 times component depth is most satisfactory 
while a width of 1.3 times element depth may be acceptable. Further 
confirmation of the observation that spacings beyond 8 mm do not sig- 
nificantly affect heat transfer is seen in Figures 64 through 68. These 
show the Nu* vs. Gr* data along with the correlation, Nu*¥ = 1.88 
(Gr*)-15, obtained in Chapter III. The correlation for an unshrouded 
surface holds very well through 14mm. At 8mm there is a slightly 
greater deviation. As anticipated, it is inapplicable for the smallest 


spacing of 6.2 mm. 


B. INDIVIDUAL ELEMENTS POWERED 

As in the absence of a shroud, individual blocks were powered to 
1.0 watt, one at a time, but with a shroud spaced 14 mm from the sub- 
strate. The resultant data is displayed in Figures 69 through 73. The 
plots follow similar trends as for the no-shroud condition (Figures 21- 
25) but show somewhat smaller temperatures. Again, the protruding 


elements have minimal effect on the unheated starting length. 
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Figure 64. Nu* vs. Gr*, With Correlation, 23 mm Channel Spacin 
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With Correlation, 6.2 mm Channel Spacin 
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Figure 68. Nu* vs 
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Figure 69. Nondimensional Temperature Excess, Front Faces 


Individually Heated Components, 14 mm Channel Spacin 
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Figure 70. Nondimensional Temperature Excess, Right Faces 


Individually Heated Components, 14 mm Channel Spacin 
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Figure 71. Nondimensional Temperature Excess, Top Faces 
Individually Heated Components, 14 mm Channel Spacin 
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Nondimensional Temperature Excess, Bottom Faces, 


Individually Heated Components, 14 mm Channel Spacin 


Figure 72. 
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Figure 73. Nondimensional Temperature Excess, Heaters 


Individually Heat 


omponents, 14 mm Channel Spacin 
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C. COMBINATIONS OF ELEMENTS HEATED 

In a manner identical to that without a shrouding wall, progres- 
sion of elements heated was made up the array while all upstream 
blocks remained powered. Additionally, block 8 was powered 
throughout and a shroud 14 mm from the substrate was in place. The 
plots (Figures 74-78) reveal that for this channel width, as with no 
shroud, the powering of downstream components has minimal effect 
on those upstream components already powered. The plots further 
show a weak effect by a heated element (block 8) on those compo- 
nents between it and the ones powered upstream, which are heated 


only by the convective flow generated upstream of them. 


D. ALL ELEMENTS POWERED, TRANSIENT 
1. Flow Visualization 

Figure 79 shows transient flow development in the x-y plane 
following 1.0 watt power step-up with a 19 mm spaced shroud in 
place. All exposures are for 20 seconds. The photographs from left to 
right are initiated at 0, 20, 40, 100, and 140 seconds following start of 
heating. In the first picture, parallel flow is developing over blocks 3 
through 8 with relatively little movement at blocks 1 and 2. By picture 
number 3, the flow follows more closely the contours of the protrud- 
ing blocks and originates further upstream. By photograph 5, the flow 
appears nearly fully developed. 

The power step-down from 0.2 watts was again for a 19 mm 


channel width and 20-second exposures (Figure 80). The pictures 
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Figure 74. Nondimensional Temperature Excess, Front Faces 


Multiple Components Heated, 14 mm Channel Spacin 


100 





SdNAL AGISHIWN VI 


B'O 





« 0 


£5 
S 
= 
S 


meee iii ia ai i i i i ia i i eS 





ee 


dhl. ae | = 


mannii BREN x 
\ LULL 9 b-- Vi 


=A, See 


NLL 88: su 











Nondimensional Temperature Excess, Right Faces 
omponents Heated, 14 mm Channel Spacin 


Multiple 


igure 75: 


FO 





SdNaAL dOL WW ¥1 _ 
ay x 
80 FO cO 





— GLH TIV 4 
(LIE Weyl io 
CUI PYS=1 x 
“CGH m H+ 
(LU QReAle vy 
GALT P9e i O 

















——-— — — -——-—- 


mm ii i ee eB eS eK KK 


G2 0SZ S2z 002 SZI 


O0£ 


eze 


6-7 «L/LLG 


Figure 76. Nondimensional Temperature Excess, Top Faces 


14 mm Channel Spacin 


Multiple Components Heated 
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Nondimensional Temperature Excess, Bottom Faces, 


Figure 77. 


Multiple Components Heated, 14 mm Channel Spacin 
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from left to right were initiated at 0, 60, 200, 320, and 460 seconds, 
respectively, following cessation of power. As with infinite channel 
width, the flow characteristic persists over the entire period observed. 
Only in the last picture does the flow begin to follow less closely the 
contour of blocks 7 and 8 and originate farther downstream. 
2. Quantitative 

The transient temperature response curves for powers of 1.0 
and 0.2 watts and a shroud spacing of 19 mm (Figures 81 and 82) fol- 
low exactly the same pattern as those for no shroud (Figures 12 and 
13). For the spacing examined, the shroud does not affect the trans- 


port evolution. 
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Transient Temperature Response 


Following Initiation of Heatin 


Figure 81. 


19 





Los 7e OF 


& 


Paps 
Ca 


C4 
‘a sewelae. 


4° oo wm 


L | 


ery 





in in in 
w @ IN ~ (0 
onal =! c=) T=! = 


16 


80 120 160 200 240 280 
TIME CSEC) 


40 


CJ> syNnLVaadWwsal 


Figure 82. Transient Temperature Response 
Following Cessation of Heating, 0.2 Watts 
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V. NCLUSION 


Flow visualization in concert with quantitative analysis provided a 
number of significant determinations: 
e The convective flow was three-dimensional. 
e The flow was laminar. 


e The thermal layer was imbedded within the momentum exchange 
region. 


e Heated elements have minimal effect on those elements, 
upstream, which are already heated. 


e There is no gain in heat transfer for channel widths greater than 
two times component depth. 


e A correlation has been found to predict heat transfer rates for the 
subject geometry. 
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VI. RECOMMENDATIONS 


It is suggested that the following areas of study be experimentally 
explored: 
e Using the same test surface: 
- explore transient response for other block faces 


- with the traversing thermocouple probe, determine variation of 
temperature with position (in the thermal layer) for all blocks 


e Construct new test surface and rerun the same battery of experi- 
ments to determine significant changes when: 


- heated elements are flush mounted 


- heated elements are spaced closer in the x-direction. 


APPENDIX A 
FTWARE 


TEMPERATURE ACQUISITION PROGRAM 


1@ ee ee ee ee ee Ds UU ab 

22 ot TEMPERATURE ACQUISITION PROGRAM it 

32 V1 (STEADY STATE) V1 

4Q i 

52 REAL Volts( 60) 

6@ REAL Temp(S3) 

72 CREATE BOAT “THES _: ,700,0,2" ,48 

8 ASSIGN @Path3 TO "THES: ,700,2,2" 

3Q PRINT " BLOCK #8" 

1O@ PRINT " " 

11 OUTPUT 709;"°CONFMEAS OCV,100,USE 0" 

12@ ENTER 709; Volts( 60) 

130 OUTPUT 709; "“CONFMEAS DCV, 100-105, USE " 

140 FOR I*@ T0 5 

150 ENTER 7O9;Volts(I) 

162 Temp( I )=.00067397+( 25825. 1328#Volts( I) )-( 607789. 2467#( Volts(I)*Volts(I1I)))-<( 
21952034. 33649(Volts(1)°3))+( 8370810996. 18749( Volts(I)*4)) 

17 PRINT “T. C. &°%3141." Volts 0.C. “;Volts(1),° Temp. DEG. € “3 Teaaqe 
180 NEXT I - 

19 PRINT " “ 

2 OO PRINT “ BLOCK #7" 

21@ PRINT “" " 

220 OUTPUT 709;"CONFMEAS DCV,106-111,USE 0" 

230 FOR I=6 TO 11 

240 ENTER 709;Volts(I) 

252 Temp( I )=.00067397+(2582S5S.1328*#Volts( I) )-( 607789 .2467#( Volts(I)*#Volts(I)))-<( 
21952034. 3364#( Volts(I)*3) )+( 8370810996. 1874#( Volts(I)*4)) 

260 PRINT “T. C. &"31+1," Volts 0.C. “s:Volts(1),” Temp. DEG. C Filename 
270 NEXT I 

282 PRINT “ " 

290 PRINT “ BLOCK #6" 

300 PRINT “ " 

312 OUTPUT 709;"°CONFMEAS OCV,112-117,USE @" 

320 FOR [212 TO 17 

330 ENTER 709;Volts(I) 

342 Temp( I) =. 0@06797+( 25825. 1 328#Volts( I) )-( 607789. 2467#( Volts(I)#*Volts(1)) )-¢ 
21952034. 3364#(Volts(1I)*3))+( 8370810996. 1874e(Volts(1)°4)) 

350 PRINT "T. CC. &"3I141," Volts O0.C. “s3Volts(I),“Temp. DEG. C “3 Temp(l) 
362 NEXT I 

370 PRINT " “ 

380 PRINT “ BLOCK #5" 

330 PRINT “ “ 

400 OUTPUT 709;"CONFMEAS OCV 118-119 ,USE Q" 

412 FOR I=18 TO 19 

420 ENTER 709;Volts(I) 


La 


4 3@ Temp( 1)=.0006797+( 25825. 13288Volts(1) )-( 607789 .2467#(Volts(I)*Volts(I)))-( 
21952034. 3364*(Volts(1)*3))+( 8370810996. 1874#( Volts(I)*4)) 

440 Pounieets ee & siti Volts D.C. “;Voltst(]),“Temp. DEG. C ";Tempt!) 
450 NEXT I 

462. OUTPUT 709;"CONFMEAS DCV,200-203,USE 9" 

47 FOR I=20 TO 23 

480 ENTER 709;Volts(1) 

490 Temp( I )=.00@6797+( 2582S. 1328#Volts( I) )-(607789.2467#(Volts(I)#Volts(1)))-¢ 
21952034. 3364#(Volts(1I)*3) )+( 8370810996. 1874#(Volts(I)*4)) 


500 PRINT “fT. C. 8° ;1+1,° Volts D.C. “;Volts(I),"Temp. DEG. C “;Temp(l) 
51 NEXT I 

22 PRINT" “ 

S530 PRINT " BLOCK #4" 

540 PRINT * * 


550 OUTPUT 709; "CONFMEAS DOCV 204-209 ,USE 0" 

560 FOR I1=24 TO 29 

S70 ENTER 709;Volts(I) 

580 Temp( I )=.0006797+( 25825. 13Z28#Volts( I) )-( 607789. 2467#(Volts(I)*#Volts(I)))-¢ 
21952034. 3364#(Volts(I)*3))+( 8370810996. 18748( Volts(I)*4)) 


S90 PRINT “T. C.. 83141," Volts O.C. "“3Volts(I),“Temp. DEG. C “s;Temp(I) 
620 NEXT I 

610 PRINT " “ 

620 PRINT “ BLOCK #3" 

630 PRINT " “™ 


640 OUTPUT 709;"CONFMEAS DCV ,210-215,USE @" 

650 FOR I[=3@ TO 35 

660 ENTER 709;Volts(I) 

679 Temp( I )=.0006797+( 2582S. 1328#Volts( I) )-(607789.2467#(Volts(I)#Volts(1I)))-¢ 
21952034. 3364#(Volts(1I)*3) )+( 8370810996. 1874#( Volts(1)*4)) 


689 PRINT "T. C. "3141," Volts O.C. “;Volts(I),"Temp. DEG. C "“;Temp( 1) 
690 NEXT I 

70@ FOR J=1 TO 14 e 

710 PRINT " " - 

720 NEXT J 

730 PRINT " BLOCK #2" 

74@ PRINT " “ 


750 OUTPUT 709;"CONFMEAS DCV ,216-219,USE 0” 

760 FOR I=36 TO 39 

770 ENTER 709;Volts(I) 

789 Temp( I )=. 0006797+( 25825. 1328#Volts (I) )-(607789.2467#( Volts(I)#Volts(I)))-<( 
21952034. 3364#(Volts(I)*3))+( 8370810996. 1874#( Volts(1)*°4)) 

79@ PRINT "T. C. #"%314+1," Volts O.C. "“sVolts(I),“Temp. DEG. C “s:Temp( I) 
802 NEXT I 

810 OUTPUT 709; "*CONFMEAS DCV , 300-301 ,USE 2" 

820 FOR 1242 TO 4! 

830 ENTER 709; Volts(I) 

B40 Temp( I )=. 0@06797+( 25825. 1328#Volts( I) )-(607789.24678#( Volts(I)#Volts(I)))-( 
21952034.3364#( Volts(I)*3) )+( 8370810996. 18748( Volts(I)*4)) 


850 PRINT "TT. C. #°31+!,”° Volts D.C. “;Volts(I),“Temp. DEG. C “sTemp(l1) 
860 NEXT I 

eo PRINT" “ 

880 PRINT “ BLOCK #1" 

890 PRINT " " 


3 O@ OUTPUT 709;"CONFMEAS DCYV ,302-307,USE 2" 
91 FOR 1=42 TO 47 
920 ENTER 709;Volts(1) 
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930 Temp( I )=. 00@6797+( 25825.1328#*Volts( I) )-(607789.24678( Volts(I)*Volts(1)))-<( 
21952034. 3364#( Volts(1I)*3))+( 8370810996. 18748( Volts(I)*4)) 

940 PRINT "T. C. 8";1+1," Volts B.C. “;Volts¢1). “Temp. DEG 9G “siemouie 
352 NEXT I 

962 | 

970 | 

980 OUTPUT 8Path3; Temp( #) 

390 | 

10@@ =| 

1010 PRINT “ “ 

1020 PRINT “ “ 

1@30@ PRINT “ BATH TEMPERATURES (TOP TO BOTTOM)" 

1@4Q@ PRINT “ " 

1050 OUTPUT 709;"CONFMEAS OCV, 317-319 ,USE 0" 

1060 FOR [257 TO S39 

1070 ENTER 7O9;Volts(I) 

1080 Temp(I)=.0006797+( 25825. 1328#Vol ts( 1) )-(607789.24674#(Volts(I)#*Volts(1)))-¢ 
21952034. 3364#( Volts(1I)*3))+( 8370810996. 1874#( Volts(I)*4)) 

109Q@ PRINT “VOLTS O.C. "“";Volts(I),"TEMP. DEG. C “;Temp( I) 

1100 NEXT I 

1118 ENO 
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DATA REDUCTION PROGRAM 


eb ie eres en Bee eee es Ge: we ie) os 


2d | M1 
3 VI DATA REDUCING PROGRAM (BASIC) l| 
40 Mt I 
50 a 


60 REAL Temp( 48) ,Tt( 48) Qcond( 8) ,Qconv( 8) ,Qflux( 8), HC8) _Nu( 8) Nustar( 8) 
72 REAL Gr( 8) ,Ttt(8) ,At(8) Otf( 8) Ote (8) ,Ot1(8) Ott’ 8) ,Otb( 8) Oth’ 8) B18) 
8 REAL Jstar(8) ,Jp(8) ,Grstar( 8) 

90 REAL Power ,Space ,Oxpg,DOxr ,Kpg,Kr , Block 

100 REAL A,G,D0is,Try,Kh2o,Beta,Kinvisc,Pr 

110 INTEGER I,J,K,L,M,N,0 

120 | 

130 Dx pg= . 006731 

140 Oxr=.@03175 

150 Kog=. 1421 

16 Kr=.0389 

17 Chpht= . 00787 

180 Chpdpthe.@@610 

19 Chpuwdth=.22 388 

2 OO Prim=(2#(Chpht+Chpwdth) )+(2«(24#( Chpht+Chpdpth) ))+(2#(2#( Chpwdth+Chpdpth) ) ) 
212 Al=Chpwdth*Chpht 

220 AZ=Chpwdth*tChpdpth 

2 30 A3=Chpht*Chpdpth 

240 Atotal=A1+(2.#AZ)+(2.#A3) 

250 Atotala=(2.#AZ)+(2.#A3) 


260 6=9.81 

270 PETEPEEEEE EEE EPEC EEE EPO EPEC EEC U EEE EEE Eee eee 
280 14 CALCULATE CHARACTERISTIC LENGTH | 
290 ee 
3@@ Lbar=Atotal/Prim 

31 ! 


320 R1=DOxpo/( Kpg#Al ) 

330 R2=Dxr/(KreAl) 

340 I 

350 PRINTER IS 1 

360 PRINT “INPUT RUN #, INTEGER ONLY" 
370 INPUT Try 

380 PRINT " “ 

390 PRINT " “ 

400 PRINT “INPUT POWER VALUE IN WATTS" 
410 INPUT Power 

A20 PRINT “ “ 

430 PRINT “ “ 

440 PRINT “INPUT BLOCK #°S POWERED" 
450 INPUT Blnp$s 

A460 PRINT " “ 

470 PRINT " “ 

480 PRINT “INPUT SPACING IN MM" 

490 INPUT Space$ 

500 PRINT " “ 


lls 


S1@ PRINT " " 
520 PRINT “INPUT AMBIENT TEMPERATURE" 
$30 INPUT Tinf 


542 PRINT i 

550 PRINT " " 

S6Q | PRINT “INPUT NEW FILENAMES" 
570 | INPUT Neufiles 

580 i PRINS 

592 Le APRN Ve 


5Od PRINT “INPUT OLD FILENAMES" 
612 INPUT Oldfile$ 


622 I 
630 | 
640 | CREATE BDAT Newfileé,128 
650 ' ASSIGN @Road TO Newfiles 
660 | 


67@ ASSIGN @Street TO Oldfiles 

68Q ENTER @Street: Temp(e¢) 

690 FOR I=0@ TO 47 

700 J=J+I 

712 IPRINT J,Temp( I) 

720 NEXT I 

730 | 

742 | 

750 PEPUPE PEEP POPPE POPP POP ieee 
762 II CALCULATE AVERAGE SURFACE TEMPERATURES It 
770 PPP OPT PEEP EPP EPP EPP ie eee 
780 | PRINT “AVG. SURFACE TEMPERATURE" 

790 J=0 

8G FOR I=1 TO 8 

812 IF I=6 THEN 

820 At(T)=(¢( Temp( J+1)#AZ)+( Temp(.J+2 )#A3)+( Temp( J+3) 0 AZ )+( Temp( J+4)#A2Z))/ 
Atotala = 

830 EVoG 

840 At( I )=(( Temp( J) e#AL)+¢ Temp( J+1)#A3)+( Temp( J+2)#AS)+( Tempo( J+ 3)# AZ )+( Te 
mp( J+4)#A2))/Atotal 

850 END IF 

862 ' PRINT At(I) 

872 J=J+6 

880 NEXT I 

890 | 

9B2@ POPC ECVPTL UOTE O OP Pe ee ii 
912 1 | CALCULATE TFILM VI 
920 PEPPP PET EOE E EEE EC EEP UOTE E 1 CLL Yee eet toe 
930 J=1 

942 Thar=(At( Ji+At( J+] tAt( J+2+At( J+ 23 +At  Jt+4) +At( J+S)+At( J+6)+At( J+7))/8.0 
950 TfRalmc=(( Tbar+Tinf )/2.0) 

960 TfRilm=Tfilmc+273.15 


972 | PRINT ~TRIEM DEGREES Ci ifilme 

38 | PRINT TELLS tfile 

99@ SU USSU SU S| 
1G@2O CALCULATE WATER PROPERTIES II 


116 


101 
102@ 
1930 
104@ 
1050 
1060 
1070 
198 
1039 
110 
111 
112 
1130 
114@ 
115@ 
1160 
1177 
118 
119@ 
12@®@ 
121 
1220 
1230 
1240 
125 
1260 
127 
128 
1290 
1300 
1310 
132 
1330 
1340 
1350 
1360 
137 
1380 
1390 
1400 
141® 
1420 
1430 
1442 
1450 
1469 
1472 
1480 
1490 
1500 


PE) UD DUS UeCIRURUIUISTISTSLISSUISUSU STP USUSSCI HG ISUIS GSS UES 5-8 FTG Tis UI TSU RSL CR 
Beta=(-7.9448S548E-8e( Tf ilm*2))+(5.7356479E-SeTfilm)- .0097810563 
Spvol#(4.69699E-9e( Tfilm*Z) )+(-2.S374S5E-6eTfilm)+. 001341903 
Dynvisc#=( 3.234851 1E-7e( Thilm*2) )4+(-.00O21474487#TF11m)+.036166792 
Kinvisc#Spvol#Dynvisc 

Kh20=(1.418181818E-3*Tfilm)+.1866 

Pr=(4.6S5706E-38(Tfilm*2) )4(-2.922094eTFf11m)+463. 3319 


Rise ces ote sual) <8) Os Le mene UM ears kelete: (is) isi cui,’ (sp lel’ fa) le: MgO!) «eee. ie) te) Gel et I) lere tele ce ese en ls 


| PRINT "TEMP-TINF” 
FOR I12@ TO 47 
IF I*S THEN 
Tt(1)=888. 888 
ELSE 
IF I=30@ THEN 
Tt(I)=888. 888 


ELSE 
Tt( 1) =Temp( 1)-Tinf 
END IF 
END IF 
| PRINT Tt¢ 1) 
NEXT I 
OOOO eee eee eee eee eee eee 
1! CALCULATE CONDUCTION LOSSES THROUGH THE TEST SURFACE li 


eo VEre Meiare he URN iaig ten ONS) e:elgltal ven Biers (@:) o> Cee (Wig os) ce. te woos, fe a, et ee le wt me) Gs fe, ee ee mm ese fal pee me cee Iw ee ey ce Se ae Te 


PRINT “QCOND" 
L=5S 
FOR I=1 TO 8 
IF I=1 THEN 
Qcond( 1)=.01S*#P ower = 
ELSE aad 
Qcond( I )=Tt( lL )#¢( 1.0/¢(R14RZ)) 
ENO IF 
| PRINT Qcond(I) 
L=L+6 
NEXT I 


oo Ola eta alan ls cua atau! Jee lkca: 


of A Oe oA nee Ms OR Ty Ue ais te ony) cat Ge fale ce ny ie ee eg wea eS Ket se te om 


| PRINT “QCONV"” 

FOR I=1 TO 8 

Qconv( I )=Power-Qcond(I ) 
Qf lux(I)=Qconv(1)/Al 
NEXT I 


iin ei eae ia vNMEn Get [i te Sia) Ae REMC ws eM) sey iets pe et ‘me eh ea bw 


| CALCULATE TEMPERATURE SCALING FACTOR 


: | : ' 8 . | . : oo. . oo. is ree Biah 19p ete: ons: eo 8 . | . . : . | | . » 8 . _ 8 « * . * . . . . 


| PRINT "Te" 
FOR I=1 TO 8 


He Oe 4 


151 
152 
153 
154@ 
155@ 
156 
1S7@ 
1S8@ 
15390 
1600 
161 
162@ 
163 
1642 
1659 
1662 
167@ 
168 
169 
17O@ 
171@ 
172@ 
173 
1742 
175@ 
1762 
177 
178@ 
179@ 
180 
181@ 
1820 
1830 
184 
1850 
1862 
187@ 
1880 
1890 
19 
191@ 
192@ 
19 3@ 
194@ 
195 
13960 
197 
198 
199@ 
2000 


Tstar(I)#Qflux(I)*tLbar/Kh20 
| PRINT TSTAR(I) 


NEXT I 


FOR [1 


TO 8 


ae re ee te ek my el ee eee mi 78 


Grstar(I )=GeBetarQflux(1)e(Lbar*4)/(KhZ0*( Kinvise*2)) 


NEXT I 


a tw al ay iw 


FOR I=! 


TO 8 


Teri Mnromrre se Ser Cor a> ois oo ri ey ie tr ok oe Ho Os 


et er Gee fe Sy Yan Yee See eet Soe te A ee SOE Yt fee tear Ws Geel itte nr Wi Be CF Oe ry a ih It skye a og 


HCD) =Of Fux 7A C1 = ant) 


NEXT I 


FOR [=] 


Oss 


CHM tee nitte tte. ete oe: 


Nu( I)*H( I) *bLbar/KhZo 


NEXT I 


FOR I=1 


Js=9- 


I 


Tens 


Sah at) Gm ee ene etc ce. pie: ee ee ian ae eee ces ee ome ee et oo 


Nustar(I)=*Nu( I )#Js*( 1/6) 


NEXT I 


Pte US A a eee 


CULV EMEA T TET a ee a 
a 


PRINTER IS 71 


{ 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


FOR I=} 


USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


“BLOCK #",“QCOND” ,“QCONV” ,"GR 8" ,“T* ">, H 7 NU ao NUe 
TO 8 


GENERATE OUTPUT DATA FILES 


— 


“""RUN NUMBER” 15x 00"; Try 

““"“POWER IN WATTS" ,11X,Z.D";Power 
“""SPACING IN MM"" 12X,10A"; Spaces 
“"" BI OCK #'s POWERED”, 8X,10A";81nps 
“““AMBIENT TEMP. DEG. C"",SX,0D.DD";Tinf 
“""E TEM TEMPERATURE" 9X ,00D0.D0";Tfilm 
"“"ETLM TEMPERATURE, C"" ,6X,DD.00";Tfilme 
“""THERMAL COND. OF WATER", 3X,Z.00D";Kh2Zo 
“EXPANSION COEFFICIENT,B"" 2X, 3DE"i Beta 
“““KINEMATIC VISCOSITY (6X. SDE Gh anmvacc 
“"““BRANDTL NUMBER" 11x ,00.00°;Pr 


UUs 


rit 


Caer oe Tne fect J 


——- 


o » @ 


* * @ 


*) ive 


2010 
2O2O 


Block=9-I 
PRINT USING “3X,0,3X,S3DE ,3X,5Z.D000, 3X ,S3DE ,2X,S3DE,3X,S3DE,2xX ,S3DE,2x,S3 


DE" :Block ,Qcond( 1) ,Qconv(1),6rstar(I),Tstar( 1), HCI) ,Nu(I) ,Nustar(1) 


2030 
2040 
2050 
2060 
2070 
2@8O 
2090 
2100 
Z11@ 
Z12@ 
2130 
2142 
2150 
2160 
217 
2180 
219@ 
2200 
221@ 
2220 
2230 
224@ 


NEXT I 

PRINT >” 
PRINT " “ 
PRINT " " 
PRINT " “ 


O00 (on eee eh ease) cecet lekce Mite ince) lis ce! a) se! oe) 6 se 


smsne eters cal | een ir sorte! ein el en ete te ise ces et ce: Se! oe) 1S) ie eelntsh i Mie ce oe ee ce ce ey Ee ow ee CH We iw ee 1x 


K=@ 
FOR N=! TO 8 
Dt f(N)=Tt(K ) 
Dtr(N)=Tt( K+) 
DtlCN)=Tt(K+2) 
Dtt(N)=Tt(K+3) 
Dtb(N)=Tt(K+4) 
Oth(Nd=Tt( K+S) 
K=K+6 
NEXT N 
PRINT BLOCK 8 . OTF ., DIR’, DIL’ , “OTT”, °OTB”, “OTH” 
FOR I=] TO 8 
Block =9-I 
PRINT USING “3X _D. 3% .S3DE.3X,S30E ,3X.S30€ , 3X,S3DE ,2X .S3DE ,2xX ,S3DE"; Block , 


Greet), Dtr<1) ,Ot1¢ 1) ,.Ott(1) Otb( 1) Oth!) 


2250 
2260 
227 
2280 
2290 
2 300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2462 
2470 
2480 
2490 
2500 


NEXT I 


re ec) Meets 


FOR I=1 TO 8 
| IF I=! THEN 
| Dthn( I )=88888 
EESE 
| IF I=6 THEN 
| Otfn( I) 288888 
Otfn( I )=Dtf(l)/Tstar (1) 
Dtrn(I)#Dtr(1)/Tstar(1) 
Otin¢ 1 )=Otl¢(1)/Tstar (1) 
Ottn(1)=Dtt(l>/Tstar(1) 
Dtbn(1)=Dtb(1)/Tstar(1) 
Othn( 1 )=Dth(1)/Tstar(1) 
IF I=! THEN 
Othn( I )=88888 
BE SE 
IF I=6 THEN 
Ot fn(1)=88888 
END IF 
END IF 


= 


—_— 


NEXT I 
PRINT " "™ 
PRINT " " 
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251 PRINT "X/L" ,“DTFN” ,“DOTRN" ,“DTLN® ,“DTTN” ,“OTBN” ,“OTHN” 
252@ Xs=.795 

2530 § Gap=.@254/ . 3048 

2540 FOR I=! TO 8 

2550 Xl=Xs-((I-1)#6ap) 

2560 PRINT USING ".DO0, 3X,S3DE ,3xX,S3DE , 3X ,S3DE , 3X, S3DE ,2X,S3DE ,2X , S3DE";X1 , Otf 
n(I) Oten(l) Otln( I) Ottn(l) ,Otbn(l) ,Othn( 1) 

257@ NEXT I 

'258@ PRINT" “ 

Z2S9® PRINT " “ 

2500 PRINT “OTR = RIGHT TT. €. = TING etc.s 

261@ PRINT “Te* IS TEMPERATURE SCALING FACTOR" 

2620 PRINT "GR # IS MODIFIED GRASHOF NUMBER" 

2630 PRINT "H IS HEAT TRANSFER COEFFICIENT” 

2640 PRINT "NU # IS NUSSELT NUMBER" 


2652 | OUTPUT @Road;Qcond( #) ,Qconv(#) ,Gr(*#) ,Ttt(#) ,Bl(#) OtF(#) Dtrie), 
Dtl(e#) Ott(#) Otb(*#) Oth #) 

2660 ! 

2670 '!' NOTE THIS PROGRAM IS USED IN CONJUNCTION WITH THE 

2680 1! THE DATA AQUISITION PROGRAM AND ACCOUNTS FOR THE 

2630 1! JNOPERATIVE THERMOCOUPLES ON BLOCKS 388 

2702 | 

Z71@ END 


120 


THERMOCOUPLE VARIATION WITH TIME PROGRAM 


1 20 DT CSET Tp CCT aCe Te SCC CPs ea a 


Base ease saael tes tal aie ee neers ete 65 UN Sher ve ee ice! et Sen ia aha 


2@ | TEMPERATURE TRANSIENT RESPONSE PROGRAM I 
3@ LUST AT A DART TS iS TS TE ek US ne a a a a i 
4@ REAL 172(575) ,173(57S) ,74(S75) , 79(975) , T6¢575) , 770575) , T8¢5795) 

5 REAL ¥2(S575) ,V3(S7S) ,V4(S7S) ,¥VS(S7S) , ¥6( S75) ,V7(S57S) ,V8( 575) 

6@ REAL 6(S75) 


7 N2S75 

82 PRINTER IS 701 

3@ CREATE BDAT “TRANS >, 708,9,1" ,S@2 
1O@ ASSIGN @Traill TO “TRANS: ,700,9,1" 
112 PCCP PET OEE EE CEP ECT PETE TEE eee 
120 BEEP 
130 FOR J=1 TO S 

14 OUTPUT 709;"CONFMEAS OCV, USE Q” 
1S5@ ENTER 709;V1(J) 

16 OUTPUT 709;"CONFMEAS OCYV,  ,USE @” 
17@ ENTER 709;V2( J) 

18 OUTPUT 709;"CONFMEAS OCV, USE @" 
19@ ENTER 709;V3( J) 

200 OUTPUT 709;"CONFMEAS OCV, USE Q” 
21 ENTER 709;V4( J) 

222 OUTPUT 709;"CONFMEAS OCV, (USE 2" 
230 ENTER 709;VS( J) 

242 OUTPUT 709;"CONFMEAS OCV, ,USE " 
252 ENTER 709;V6&( J) 

260 OUTPUT 709;"CONFMEAS DCV,  _,USE @" 
27 ENTER 709;V7( J) a 
280 OUTPUT 709;"CONFMEAS OCV, (USE @” 
290 ENTER 709;V8& J) 

30@ NEXT J 

310 BEEP 

320 FOR J=6 TO N 

330 OUTPUT 709;"CONFMEAS OCV, (USE Q" 
34@ ENTER 709;V1(J) 

350 OUTPUT 709;"CONFMEAS OCV, (USE @" 
360 ENTER 709;V2(J) 

37@ OUTPUT 709;"CONFMEAS DCV, USE @" 
380 ENTER 709;V3( J) 

39@ OUTPUT 709;"°CONFMEAS OCV,  ,USE @Q" 
400 ENTER 709;V4( J) 

412 OUTPUT 709;"CONFMEAS OCV,  _|USE Q" 
420 ENTER 709;VS( J) 

430 OUTPUT 709;"CONFMEAS OCV, _,USE 9" 
44@ ENTER 709;V6(J) 

452 OUTPUT 709;"CONFMEAS OCV, USE 0" 
460 ENTER 709;V7( J) 

472 OUTPUT 709;"“CONFMEAS DCV, USE 0" 
480 ENTER 709;V8( J) 

49@ IF J=N THEN 

SOO BEEP 2450.97#Tone,!.S 

S10 END IF 

S2Q NEXT J 


Ze 


530 
540 
550 
560 
57@ 
580 
59@ 
bO2 
610 
629 
630 
640 
65 
660 
670 
68 
690 
70@ 
71 
72 
730 
74 
750 
76@ 
77 
780 


A= .0006797 
B=25825. 1328 
C=-607789.2467 
D=-21952034. 3364 
E=8370810996. 1874 

FOR J=! TON 

IF J<S THEN 

G( J) =@ 

ELSE 

G(J)=(J-S)*Tim 

1! Tim IS TIME BASED ON TIME PROGRAM RUNS FROM 2ND BEEP TO TONE DIVDED 

' BY N-S, E.6. FOR 7 THERMOCOUPLES Tim IS .74S, FOR 8 IT IS .852 
END IF 
NEXT J 
PRENT =. 
PRINT ° " 
FOR J=1 TON 
TiC J) #At(BeVid 5) 40°C #CVid 5) 72) 4008 V1 S)739) 4 Ee CV J) *4)) 

T2¢ J) =A4(B0V205))4¢C8(V2059°2)04(¢080V205)°35) 408 080V¥205)°4)) 
T3( J) 2A4( BeV3( 5) )4°¢C8¢V30I)°2) )4¢D8(V305)°3) 4(E0(V50J)*4)) 
T4( J) =A4+( BeV4ad 5) )4+¢C#(V4(05)°2) 4¢08(V4(0J)°53))408E 0 (V4 J)*4)) 
TS( J) =A4( BeVS(0 5) 4¢C#(V5059)°2) )4¢00°0V505)°3) 40F 00 V50I)°4)) 
T6( J) =A+(BeVG( J) )4¢C#(VG(5)°2) 14¢ D8 VEC J)°3) 406 *#(VE0I)74) ) 
T7( J) =A+(Bev7( J) )4+¢C8CV7059°2) )4¢ 000470 5)°3) 4CE eC V705)*%4)) 
T8( J) 2A+( B#V8( J) )4¢C#(V8¢I)72))4¢08(V8(5)°739) +0 Ee VBC I)%*4)) 

PRINT USING “1X,000,!1X,40E,!1X,40E,1xX,4D0E,1xX,40E,1X,40E,1xX,4DE ,1X ,4DE,1x,4 


DE“ 3G¢ 3). T1(J).12¢3) 1300) , 1407). TSC). GC) bee 
799 NEXT J 
82 OUTPUT @Trail1;G(*) TI #) T2(#) TH) Ta) TSC *) TEC#) T7(*) TECe) 


819 


END 


PAL 


APPENDIX B 
EQUATIONS FOR DETERMINING FLUID PROPERTIES 


2 
B = -7.945E-10T-,,, + 5.736E-5 T py — 9-781E-3 


v =4.697E-9T_,, — 2.537E-6T,,,,, + 1.342E-3 


= 3.235E-7 T_,, — 2.147 E-3 T,,,, + 0.03617 


k = 1.42E-37,,,, + 0.187 


Pr =4.66E-3T , 7, — 2-92 T poy, + 463 


Curve fits were determined for range of temperatures from 280 to 


300 K. 


123 


APPENDIX C 
AMPL ALCULA 


The calculations were based on temperatures derived from a run 
with all blocks powered to 1.0 watt, no shroud in place. When 


appropriate, values are calculated for block number 6. 


A. CHARACTERISTIC DIMENSIONS 


2(0.00787 + 0.02388) + (2( 2(0.00787 + 0.00610))) 


= 0,2 
+ (2 2(0. 02388 + 0.00610))) 0.2393m 


Perimeter = 


A, = (0.02388)( 0.00787) = L879 x 10° m? 
A, = (0.02388)( 0.00610) = 1457 x 10m? 
A, = (0.00787)( 0.00610) = 4.80 x 107° m? 


A... = 4.879 x 10°) + (2)(1.457 x 107“) + (2)(4.80 x 10°) = 5.75 x 1009s 


TOTAL 


A =(5.75x 10° — 1879 x 105.) = siege loners 


TOTALA 


- §75x10— 


6) 
jo 0.2393 =2.40x 10 m 
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B. SURFACE AVERAGED TEMPERATURES 


(19.82)(1.879 x 107“) + (19.48)( 4. 80 x 107°) 
+ (19.49)( 4, 80 x 10°) + (20.62)(1.457 x 107°) 
+ (20.26)( 1 457 x 107°) 


= 20.08 °C 
(5.75 x 107°) 


T ave (8) = 


(19.41)(1.879 x 107“) + (18.99)( 4, 80 x 107°) 
+ (19.05)( 4.80 x 10°) + (19.88)( 1457 x 107’) 
+ (20.08)( 1 457 x 10°*) 


T 
(5.75 x 107’) 


AVG (7) = 


= 19.64°C 


(19.16)( 1.879 x 10°“) + (18.67)( 4.80 x 10°) 
+ (18.61)( 4,80 10°) + (20.00)( 1457 x 107‘) 
+ (19.57)(1.457 x 10°) 


T 
(5.75 x 107’) 


AVG (6) — 


= 19.40 °C 


(18.96)(1.879 x 107°) + (18.33)( 4, 80 x 107°) 
+ (18.47)( 4, 80 x 10°) + (19.68) (1.457 x 107’) 
+ (19.37)( 1.457 x 107’) 


T 
(Sas 10r ) 


AVG (5) ~ 


= 19.16°C 


(18.65)( 1.879 x 107“) + (18.30)( 4, 80 x 107°) 
+ (18.39)( 4, 80 x 107°) + (19.75)( L457 x 107°) 
+(19.06)( 1.457 x 10°’) 


T 
(5.75 x 10°’) 


AVG (4) — 


= 18.99 °C 


129 


hilton 10°°) + (18.10)( 4. 80 x 107°) | 
+ (19.18)(1.457 x 10°) + (18.87)(1.457 x 107*) ; 
fuses a ee 18.77°C 


(18.07)(1.879 x 107“) + (17.40)( 4. 80 x 10°) 
+ (17.79)( 4.80 x 107°) + (18.85)(1.457 x 107’) 
+ (18.56)( 1457 x 107’) 


lave (2) (5 oe 10") = 13.32% 
(17.73)( 1.879 x 107°) + (17.25)( 4. 80 x 107°) 
+ (17.25)( 4,80 x 107°) + (18.60)(1.457 x 107°) 
+ (17.78)( 1.457 x 10°*) ; 
owen = (5 95 5 10") = 17.89 @ 
C. TFILM 


Ff = 20. 08 + 19.64 + 19.40 + 10+ 18.99 16 eae 19. 03°C 


1G = 16.13°C 


FILM 


_ 19.03 + 13.23 
Bs Z 


D. TEMPERATURE EXCESSES 


Using data from block 6, temperature excess is evaluated as: 


AT, = 19.16— 13.23 oon 


AT, = 18.67 — 13.23 = 5.44°C 
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AT, = 18.61— 13.23 = 538°C 
AT, = 20.00 - 13.23 = 677°C 
AT, = 19.57 — 13.23 = 6.34°C 


AT,, = 24.20 - 13.23 = 10.97°C 


CONVECTED HEAT FLUX 


Ta 0. 984 — = 5240” - 
L879 x 10 m 


TEMPERATURE SCALING FACTOR 


(5240)(2.40 x 107°) 
aa, i et OSS a @) 
om (0,597) ere 


MODIFIED GRASHOF NUMBER 


4 


Gp += (9-8D(163 x 10°°)(5240)(2.40 x 107°) 


Seni 
(0.597)(112 x 107°) 


HEAT TRANSFER COEFFICIENT 


(5240) 


h= (79.40 — 13.23) 


w 
= §49——_ 
eK 


127 


I. NUSSELT NUMBER 


 (849)(2.40 x 107?) 


NS (0.597) a 


E26 
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